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^ ■ Abstract 

, There are two competing descriptions of nematic liquid crystal dynamics: the 

pP | Ericksen-Leslie director theory and the Eringen micropolar approach. Up to this 

day, these two descriptions have remained distinct in spite of several attempts 
to show that the micropolar theory comprises the director theory. In this paper 
we show that this is the case by using symmetry reduction techniques and in- 
troducing a new system that is equivalent to the Ericksen-Leslie equations and 
includes disclination dynamics. The resulting equations of motion are verified to 
^ | be completely equivalent, although one of the two different reductions offers the 

possibility of accounting for orientational defects. After applying these two ap- 
proaches to the ordered micropolar theory of Lhuiller and Rey, all the results are 
i \ eventually extended to flowing complex fluids, such as nematic liquid crystals. 
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1 Introduction 

The Ericksen-Leslie equations for the dynamics of nematic liquid crystals are widely ac- 
cepted and experimentally well-characterized through various measurements [deGennes(1971), 
Chandrasekar(1992), Rey & Denn(2002)]. However, when orientational defects (discli- 
nations) are present in the system, this model fails to provide a reliable description. For 
example, in the presence of defects, the liquid crystal molecules may undergo phase tran- 
sitions, e.g., from uniaxial to biaxial, and the director field n is no longer an appropriate 
order parameter variable. 

Among the various descriptions that incorporate defect dynamics, the micropolar 
theory developed by Eringen [Eringen(1997)] appears promising to describe the motion 
of microfluids, including liquid crystals. Indeed, besides incorporating molecular shape 
changes into a microinertia tensor j, disclination dynamics is encoded in the so called 
wryness tensor 7, which is expressed in terms of (Vn) x n when defects are absent 
[Eringen(1997)]. 

However, while nematic liquid crystals are well known to be a typical example of 
microfluids, in spite of several attempts, it is not known how the Ericksen Leslie (EL) 
description arises from Eringen's micropolar theory. For example, the relation 7 = 
(Vn) x n proposed by Eringen [Eringen(1997), formula (11.2)] fails to return the correct 
EL equations [Leslie(1979)] as shown in [Gay-Balmaz & Ratiu(2009), Theorem 8.11] by 
two different methods (symmetry considerations and a direct computation). Thus it is 
not completely clear how 7 may be expressed in terms of the director n. 

More recent developments in the understanding of defect dynamics are provided by 
the use of reduction theory [Holm(2002), Gay-Balmaz & Ratiu(2009)], which is behind 
the gauge-theory approach [Dzyaloshinskii & Volovick(1980)]. It applies to very general 
systems since it incorporates defect dynamics in different contexts, such as frustrated 
spin glasses [Holm & Kupershmidt(1988), Dzyaloshinskii & Volovick(1980)], for exam- 
ple. In this setting, one is naturally led to consider the wryness tensor 7 as the magnetic 
potential of a Yang-Mills field (or, equivalently, a connection one-form) taking values 
in the Lie algebra so (3) of antisymmetric 3x3 matrices (usually identified with vectors 
in M. 3 ) of the rotation group SO (3). The quantity 7 is also known as 'spatial rota- 
tional strain' [Holm(2002)] and it expresses the amount by which a specified director 
field rotates under an infinitesimal displacement. Due to its tensorial nature, the gauge 
potential 7 may be conveniently expressed in terms of an appropriate basis as 



7 = 7; dx l = 7? e a dx l 
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where {e a } is a fixed basis of M. 3 ~ so (3). Then, its corresponding magnetic vector field 
is given componentwise by 

B i = e i ^{d jlk + lj x lk ), (1.1) 

where we sum over repeated indices and we have used the equivalence between two- 
forms and vector fields on physical space (see §3.4 for the coordinate-free definition). In 
the gauge-theory approach developed in [Dzyaloshinskii & Volovick(1980)], the absence 
of disclinations is given by a vanishing magnetic field B, rather than by a vanishing 
potential 7. Thus, the presence of 7 in a mathematical model must be compatible with 
EL dynamics, as long as B = 0. In the context of reduction theory, one recognizes that 
a vanishing magnetic field B = simply amounts to the homogeneous initial condition 
7 = [Dzyaloshinskii & Volovick(1980)]. If the latter condition is not satisfied, then 
the gauge-theory model would extend the EL formulation to incorporate non-trivial 
disclination dynamics. 

On the other hand, Eringen's micropolar theory does not seem to possess a gauge- 
theory formulation, since the wryness tensor (Vn) x n, as defined by Eringen, does not 
transform as a magnetic potential under gauge transformations; see [Gay-Balmaz & Ratiu(2009), 
Lemma 8.10]. Nevertheless, Eringen's theory still shares many analogies with gauge- 
theory models and the coexistence of the wryness and microinertia tensors in the dy- 
namics provides an interesting opportunity to account for the shape evolution of the 
molecules interacting with disclination lines. 

The considerations above represent the main motivation for the present work, which 
uses Euler-Poincare variational methods to provide a unifying framework for incorporat- 
ing defect dynamics in continuum systems with broken internal symmetry (e.g., liquid 
crystals) and shows that Eringen's micropolar theory comprises Ericksen-Leslie dynam- 
ics. This is done upon noticing that taking the gradient of the relation 

n(x,t) = x(x,t)z, 

relating director dynamics to the dynamics of the rotation matrix x( x )^) £ 5*0(3) in 
EL theory, immediately leads to 

Vn = (V%)z = (V%)x _1 n. 

Here z := (0, 0, 1). Then, one observes that the new variable 

1 = -iyx)x~ l 

is precisely a connection one form taking values in so(3) [Holm(2002), Gay-Balmaz & Ratiu(2009)]. 
It is straightforward to see that analogous relations hold independently of the order pa- 
rameter space. Then, upon using the isomorphism so (3) ~ M 3 given by a k = —e k ^ l aji, 
one can simply replace the relation Vn = n x 7 into the EL equations to account for the 
potential 7 as an extra dynamical variable. Notice that, although the latter relation is 
satisfied by the choice 7 = (Vn) x n, this expression is only defined up to a component 
parallel to n. Thus, 7 cannot be entirely expressed in terms of the director n and it 
needs to be specified by all three columns of the matrix x(x, t). 
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The second key observation is that a different symmetry reduction of the same mate- 
rial Lagrangian yields a new set of equations for nematodynamics. We show that these 
are completely equivalent to the original Ericksen-Leslie equations. However, this new 
system allows for the description of disclinations, something that the Ericksen-Leslie 
equations could not handle, as discussed above. 

As we shall see, all the above considerations hold regardless of the background fluid 
motion and they are a peculiar feature of the micro-order. Thus, we shall mainly confine 
our treatment to motion-less liquid crystal continua in order to emphasize the high points 
of the discussion. The extension to flowing fluid systems will be presented briefly at the 
end of this paper. 

Plan of the paper. This paper starts (Section 2) by showing how reduction theory 
can be applied to Ericksen-Leslie nematodynamics in two different fashions, thereby 
producing two different sets of equations of motion. The resulting dynamical systems are, 
however, completely equivalent. In Section 3, these two equivalent reduction methods are 
then formulated in a general context, for an arbitrary order parameter space. Momentum 
map properties are presented in detail for the two constructions, which are then specified 
to micropolar continua. In Section 4, Eringen's theory of micropolar media is shown 
to comprise Ericksen-Leslie nematodynamics. This requires a specified choice of the 
micropolar free energy, which in turn reduces to the Frank energy under the assumption 
of uniaxial molecules. While Section 5 deals with the Lhuiller-Rey theory of ordered 
micropolar continua [Lhuiller & Rey(2004)], Section 6 extends all the results to liquid 
crystal flows, thereby showing how the hydrodynamic Ericksen-Leslie equations possess 
a micropolar formulation. 

2 Two equivalent reductions for nematic systems 

This section develops the guiding example of this paper, i.e., the dynamics of nematic 
media. In particular, this section shows how the reduction producing EL nematodynam- 
ics is accompanied by an equivalent reduction procedure that naturally incorporates the 
connection 7 = — (Vx)x" 1 as an extra dynamical variable. The latter construction will 
be presented after the following review of the reduction underlying EL dynamics. 

2.1 Reduction for the Ericksen-Leslie equations 

The reduction process producing the EL equations has been widely explained in [Holm(2002), 

Gay-Balmaz & Ratiu(2009)]. This process starts by identifying the configuration space 

of a nematic continuum with the space J-'iV, SO (3)) of 5"0(3)-valued scalar functions on 

the domain T> C M. 3 . Then, one makes use of the Lagrangian [Holm(2002), Gay-Balmaz & Ratiu(2009)] 
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where, usually, n = z (although it can be an arbitrary director field n (x)) and J is 
the micro inertia constant. Here the free energy is given by the Frank expression 

F(n, Vn) = K 2 (n • curl n) +7^11 (div n) 2 +^K 22 (n • curl n) 2 

chirality splay twist 

+ i^33.||nx curln|| 2 , (2.2) 



bend 



(here, K 2 7^ for cholesterics and K 2 = for nematics). The free energy can also contain 
additional terms due to external electromagnetic fields. 

In order to apply the Euler-Poincare theory for systems with broken symmetry, 
see [Holm, Marsden & Ratiu(1998), Gay-Balmaz & Tronci(2010)], we write C(x,x) = 
L(x, x,n ), where L : TJ r (V,SO(3)) x J^iV, S 2 ) — > K and observe that L is invariant 
under the right action 

(X,n ) 1 y (x^i^no) 

of ip G J-(T>, SO(3)). This invariance property yields the reduced Euler-Poincare La- 
grangian 

4(f,n) = £j / / ^(n,Vn)/i, 

where P = XX -1 an d n = x n 0! where u G M 3 H- P G so (3) denotes the usual Lie 
algebra isomorphism defined by u ab = —e a b c u c . We thus obtain the following equations 

[Holm(2002), Gay-Balmaz & Ratiu(2009)] 

d 8h Six 6h 
-— = v x hnx — 

dttfi/ 5u 5n (2.3) 

9 t n + n x v = 

by applying the usual Euler-Poincare variational principle 

5 [ £i(u,n)dt = 0, (2.4) 



to 



subject to the variations Su = 77 + v x 77 and 5n = 77 x n for arbitrary 77 G so (3) 
satisfying 77 (to) — ^(ti) = 0. More explicitly, upon denoting h = — S£i/5n, one has 

JdfU = h x n 

(2.5) 

d t n + n x 1/ = 0, 

which produce the EL equations of nematodynamics [Leslie(1979)] (upon setting i> ■ n = 0) 

_ 2 ( n • h + J n ■ n + h = . (2.6) 



dt 2 v 

Notice that the projection Ju ■ n of the angular momentum Jv on the director n 
is a constant of motion of EL dynamics that can be set equal to zero because of the 
rod-like nature of the uniaxial molecules, as explained in [Gay-Balmaz & Tronci(2010)]. 
However, when disclination lines are present, molecules may change their shape (e.g., 
from uniaxial to biaxial) and the projection Jv ■ n cannot be a constant. The next 
section considers an alternative reduction for nematic systems, which can be naturally 
extended to account for disclination dynamics. 
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2.2 Alternative reduction for nematic systems 

While the EL equations are well known and widely accepted as a reliable model, this 
section presents an alternative set of equations that are completely equivalent to EL 
dynamics. These new equations have the advantage that they can be easily extended to 
consider non-trivial disclination effects, as we shall see later. 

The starting point is the same unreduced Lagrangian (2.1) producing EL dynamics. 
We notice that, if n is an arbitrary constant director, the Lagrangian £ in (2.1) possesses 
the alternative invariance property 

£(x> X) = L (xx~\ xn , -(Vx) x -1 ) 

where V denotes the usual differentiation operator and 7 := — (Vx) X -1 e 0, 1 (V,so(3)) 
is a new dynamical variable. More precisely, if n is a constant vector field, then we can 
rewrite C as 

C{x,X) = \j I 11*11 V- / F(xn ,((Vx)x- 1 )(xn ))// (2.7) 
* Jv Jv 

and consider £ as coming from a Lagrangian L = L(x, x, n 0; 7o) defined on TJ-"(D, SO (3)) x 
J-"(P,]R 3 ) x Q 1 (T>,so(3)), which is invariant under the right action 

(X> n o, 7o) (x^> ^ _1 n , ^ _1 7o^ + ^A" 1 W>) • 
In the present case, the initial value of 7 is zero, that is, we have 70 = so that 

£(X>X) = L (x,X,no,0) ■ 
Then, the reduced Lagrangian corresponding to (2.7) takes the form 

L(xx~\ Xn , -(Vx) X" 1 ) = -J I W? V ~ I F (™, -7 * n> = i 2 (u, n, 7) , (2.8) 

^ Jv Jv 

where we allow for 7 = — (Vx) X -1 £ Q 1 (T>,so(3)) to be an extra dynamical variable, 
we denote by 7 = (7i>72>73) the corresponding M 3 -valued one-form, 7^ G r2 1 ('D), 
z = 1, 2, 3, and 7 x n e ^(D, M 3 ) is defined by 

( 7 xn)«) = 7 «) xn, v x eT x V, (2.9) 

or, in local coordinates, 7 x n = ( 7i x n) dx l . In Appendix A it is shown how the Frank 
energy is written in terms of n and 7. 

It is important to notice that the expression for L = L(x, x, n o> 7o) may not be defined 
when 70 7^ 0. In this case £2 is only defined on the orbit of 70 = 0, that is, on 7 of the 
form 7 = — (Vx)x 1 - However, this does not affect the reduction process, as long as the 
expression L(x, X? n o? 0) is invariant under the isotropy group of 70 = 0. It is interesting 
to observe that this construction is identical to the reduction process occurring for the 
dynamics of polymer chains, see [Ellis et al.(2010), Gay-Balmaz, Holm & Ratiu(2011)] 
to which we also refer for more details about the reduction processes when 70 = 0. 
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Notice that, in the context of symmetry breaking [Gay-Balmaz & Tronci(2010)], the 
above reduction is no longer performed with respect to the isotropy subgroup of z, i.e., 
F{V, 5 1 ) = T(T>, SO(3))%, as it happens for Ericksen-Leslie dynamics. Rather, since 
the isotropy subgroup of 70 = is given by 5*0(3) C F(V, SO (3)), the entire reduction 
process is with respect to the isotropy subgroup 

T(V, 5O(3)) (2 , 0) = T(V, 5 1 ) n 50(3) = 5 1 

At this point, the resulting associated Euler-Poincare variational principle is 

b [ 1 £ 2 (i/,n,7)dt = 0, (2.10) 

subject to the variations bv = rj+vxr) and S(n, 7) = (77 X n, — V 1 r]) for 17 (£ ) = rf(ti) = 
0. Here f) = (Sx)x~ 1 an d V 7 denotes the covariant differentiation V 7 a := Va + 7 x a. 
Then, the affine Euler-Poincare equations are 

d bi 2 bi 2 S£ 2 ^/ b£ 2 \ bi 2 

— — = v x — + div -=- + Tr 7 x — + n x — 
at bv bv 07 \ 07 J bn 

a n (2.11) 

o t n + n x v = v ' 

. dtl + 7 x v + Vi/ = 0, 7 = 0. 

When one allows for 7 ^ 0, the reduced equations in (2.11) still make sense, thereby 
extending EL dynamics to account for disclination dynamics. Notice that in this case, 
equations (2.11) still preserve the relation Vn — n x 7 = 0, since 

d \ 

— j/x (Vn — n x 7) = . 

dt J y " 

Thus, the initial conditions 7 and no may be strictly related through the relation 
Vn = n x 7 . It is important to emphasize that the projection n • 7 gives zero 
contribution to the gradient Vno. Then, Eringen's expression of the wryness tensor 
7 = Vn x n (that is n • 7 = 0) becomes a convenient initial condition, which is not 
preserved in time. 

Thus, if we suppose that the reference director field n is constant, then (2.3) and 
(2.11) are equivalent since they are induced by the same Euler-Lagrange equations for 
£{Xi X) 011 TJ^lV, 5*0(3)). We shall verify this fact explicitly in the next subsecion. 

2.3 Compatibility 

Upon choosing n= p and 7 = — (Vx)x _1 , the induced variational principles (2.4) 
and (2.10) must be the same: 
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rt\ rt\ rti pti 

/ £ 1 (y,n)dt = S L no (x, x) ^ = b L (n0i0) (x, x) dt = b / 

to J to J to J to 
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for any variation of x vanishing at the endpoints and o"(n, 7) = (rj x n, — Vr)), where 
rj = {5x)x l - Thus, one has 



to 



so that 



where we ignore variations in since they give equal contributions which cancel each 
other because 5£i/5v = btijbv. In conclusion, isolating 77 yields 

Mi M... „ / <%\ 

n x — = n x h div h Tr 7 x — (2.12) 

on on 07 \ 07 / 

which can be used to show that (2.11) is compatible with (2.3). Indeed, we check that 
equations (2.11) still produce 

Aa. n L4(,. n) ,0, (2.13) 

since by (2.12) all terms involving Si^/S^ are orthogonal to n and give zero contribution. 
Notice that the above constant is actually a momentum map arising from the invariance 
of EL dynamics under the isotropy group TiT*, 5 1 ) of z; see [Gay-Balmaz & Tronci(2010)] 
and also (3.4) below. The conservation of n • 5£2/Su is then to be considered a conser- 
vation law arising from Noether's theorem, which is inherited from the system (2.3). 

Notice that, if 70 7^ 0, the uniaxial property u ■ n = is not preserved, since (2.12) 
is no longer true. This indicates that changes occur in the molecular configuration of 
the system. Thus, the director parameter must be replaced by a suitable inertia tensor, 
which becomes the new order parameter field. This is precisely what happens in the 
Landau-deGennes dynamics of the alignment tensor [dcGennes(1971)]. This treatment 
is the basis of the micropolar theory of liquid crystals, which was pioneered by Eringen 
[Eringen(1997)]. However, before approaching this problem, we shall show how the 
two constructions presented in this section are actually special cases of two reduction 
processes that can be carried out for any continuous medium with broken symmetry. 



3 Reductions for continua with broken symmetry 

This section generalizes the two constructions previously applied to nematic liquid crys- 
tals to arbitrary continuum systems with broken symmetry. More precisely, the rotation 
group 50(3) for the orientational order is replaced by an arbitrary Lie group O acting 
transitively on an order parameter manifold M and so (3) = M 3 by the Lie algebra of 
O. Then, if no G M is a given order parameter variable, it follows that M is the coset 
manifold O/V, where V '■— O no C O is the isotropy subgroup fixing no- This is precisely 
the same setting as in [Gay-Balmaz & Tronci(2010)]. 

For continuous media, one replaces O and M by T{V,0) and ^{D^M), where T> 
is the spatial domain of the medium, so that T(T>, O) no = T(T>,V^) C T(T>, O), where 



S 



n e M is identified with a constant function on T>. Note that here we consider the 
action of T{V, O) on J-(T>, M) naturally inducedAt the fundamental unreduced level, 
one starts with a Lagrangian functional C : TTiT), O) —> R, which is typically of the 
type 

£(x,x) = / ^(XX,X™o, V(xn ))/i 
Jo 

so that the fixed order parameter no appears in the Lagrangian density Jzf through both 
Xno and its differential V(x^o)- In all cases under consideration, the Lagrangian C 
possesses the following invariance properties: 

£(x,x) = / -3f(x,X,X™o, V(xn ))u = / ^(xx~\x™o, V(xn ))/i = M^™) (3.1) 
and, if n is constant in space, 

£(x,x) = / ^(x,X,Xno, V(xn ))^ = / ^(xx~\ X%, ((Vx)x _1 )m (x™o)) A* 

=:£ 2 (v,n, 7 ), (3.2) 

for 1/ := XX" 1 , n := X™o and 7 := -(Vx)x" 1 - Here (Vx)x" 1 )m(x^o) : TV -»■ TM is 
defined by 

((Vx)x _1 )Af(x«o)K) = (VxK)x(^) _1 ) M (x(^)no(x)) e r n(x) M, « x e t x v, 

where £m G X(M) denotes the infinitesimal generator associated to the Lie algebra 
element (Go. Note that in the formula above we have £ = Vx(v x )x( x )~ 1 £ °- Each of 
these invariance properties involves a distinct reduction procedure that, in turn, produces 
different Euler-Poincare equations of motion. As is explained below, these two reduced 
systems are compatible since they arise from the same unreduced Lagrangian £(x, x)- 

3.1 First reduction 

This reduction procedure is based on the invariance property (3.1) and it follows pre- 
cisely the same steps as in Section 2.1 of [Gay-Balmaz & Tronci(2010)] (see theorem 2.1 
therein). In particular, this reduction is performed with respect to the isotropy subgroup 
F(V, O) no = F(p, ^n )' because of the diffeomorphism 

TF{V, 0)/F{V, Oj — ► F{V, 0) x F(V, M) 
[(x, x)] ^ (xx _1 , x^o) ; 

see [Gay-Balmaz & Tronci(2010), Remark 2.5 and equation (3.1)]. 
Since the invariance property (3.1) implies 

5 ^ £( X ,x)dt = 5 I" e 1 (u,n)dt = 0, 

Jt Jt 

the Euler-Poincare variational principle for l\ involves the variations 

=^(xx~ 1 ) = v + [v,n\ 

Sn =S(xn ) = r] M (n), 
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where rj = {Sx)x~ 1 £ •f'iPi °) an d the dot notation stands for partial time derivative. 
Here the index M on 77m(^) for V G FiP-, o) and n G J r (^ > , M), denotes the infinitesimal 
generator of the Tp, 0)-action on Tp, M), which is formally given at i G D by 
v( x )M{n(x)). Thus, the resulting equations of motion read 

est! sii fsiA 

+ ad 1/ — = J — , n = u M {n), (3.3) 



where J : T*Tp, M) — > Tp, o)* is the momentum map of the cotangent lifted action 
of Tp, O) on Tp, M), which is given by [Marsden & Ratiu(1994), Theorem 12.1.4] 

(JKU) = (a n ,£ M {n)} V« n G T*Tp, M). 

Notice that Tip, C no ) -invariance of C(x, x) — n) yields, by Noether's theorem, 
the following momentum map conservation (see [Gay-Balmaz & Tronci(2010), §4]) 

d WAd^YUo, (3.4) 



dt V. V x &v t 

where i* is the dual of the Lie algebra inclusion i : T(T>, o no ) — > T(T>, o). This is the 
immediate generalization of the relation (2.13) for uniaxial nematics, for which O no = S 1 
and i(r) = (0, 0, r). The above conserved quantity is readily seen to arise as a momentum 
map J : T*T(V, O) -> T(T>, o)* o by the following computation 

(JMX) = (a x ,Co(x)) = (a x ,X«(C)} = (x" 1 "*,^)) = <^(x _1 « x ) , C) 

where x~ x a x = Ad* (a x x~ r ) and £ G ^("P, o„ ) is arbitrary. The index O on (o denotes 
the infinitesimal generator of the right Tip, O no ) -action on Tp, O). This Lie algebra 
action is given by % \-> xKO- 

3.2 Second reduction 

In this section we restrict all considerations to a given initial condition n G M C 
Tp, M) (i.e., n spatially constant) in order to perform the reduction arising from the 
invariance property (3.2). The construction in this section is based on the treatment in 
[Gay-Balmaz & Ratiu(2009)], involving affine actions of the gauge group Tp, O). 
Property (3.2) arises mainly from the following observation: 

Vn = V(x^o) = (Vx)n = ((Vx)x _1 ) M X^o = ({^x)x~ X ) M n =: ~7mW , (3.5) 

which defines the connection one form 7 := — (Vx)x^ 1 Q l p, 0). Then, it becomes 
natural to incorporate the latter quantity in the equations of motion as an extra dynam- 
ical variable. This step requires precisely the reduction given by the invariance (3.2), 
and hence we conclude 

£(x,x) = / -^(xX~\x™o, ((Vx)x _1 )m(x«o))i" = / S?(v,n,-7 M (n))ii 
Jv Jv 

=a 2 (u,n, 1 ). (3.6) 
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In this case, the reduction proceeds with respect to the isotropy subgroup of (n (x), 70 0*0) 
(n , 0), which is necessarily a subgroup of TiT), 0) n . More precisely, since J-'iV, O) no = 
F(V, Oj and F{V, O) lo=0 = O one has 

nv,o) {nofl) = T{v,oJno = o no . 

Thus, the second invariance property (3.2) leads to a reduction involving the isotropy 
group O no , which is much smaller than the isotropy TiD, O no j used in the first reduction 
presented in §3.1 arising from the invariance property (3.1). 

Notice that, upon considering the gauge action of TiT), O) on f2 1 (P, o) 

7o H- VSoV' + W , (3.7) 
the invariance property (3.2) takes the form 

£(X,X)= / ^(xX~\x™o, -(x~ 1 7o)m(x^o)) H = / n, -7w(n)) n =: l 2 {v, n, 7) 



where 70 = is a fixed initial condition. At this point, since the invariance property 
(3.2) implies 

5 I' £(x,x)dt = 5 f 1 £ 2 (v,n n )dt = 0, 



J to J to 

the Euler-Poincare variational principle for £ 2 involves the variations 

8v =s(xx~ 1 ) =v + 

5n =5(xn ) = r) M {n) 

h =^(x" 1 7o) = -»7ni(7o) = -V 7 r/, 

where V 7 A := VA + [7, A] is the covariant differential and the index Q l on t^qi denotes 
the infinitesimal generator of the affine action (3.7). 
The resulting equations of motion are 

where 

div 7 ( — — I : = div —A- — Tr ( ad*, 



£7 J ^7 \ 7 ^7 

Here J : T*T(T>, M) — > TiT), 0)* is the same momentum map as in §3.1, while K(7, w) : = 
div 7 w is a momentum map K : T*r2 1 (D, 0) — > T{T> : 0)* arising from the lifted action of 
F(V,0) on tt l (V,o). 

3.3 Compatibility of the two approaches 

Since the two approaches arise from the same unreduced Lagrangian, they are compat- 
ible. This compatibility is reflected in the following relations 

dS£ 1 ^,,6£ 1 T /^i\ 8 5£o 5£ 2 (5h\ 7 {5£ 2 \ 
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which arise from the variational principles 

5 £i(z/,n)dt = 5 / 4(^,n,7)dt = 0. 

«/ to ^ to 

In the particular case when 8i\/8v = 8I2I8V we obtain 

which generalizes the analogous relation (2.12) previously found for nematodynamics. 

So far, we considered the case in which the parameter no appears in the Lagrangian 
density Jzf only through the term xn and its gradient V(x^o)- Then we showed how 
such a Lagrangian possesses the two invariance properties (3.1) and (3.2). However, one 
can consider the more general case of an invariant Lagrangian of the type 



£(x,x) = / =^(x,X,X^o, Vx)" = / =^(XX ,X™o, (Vx)x ) A* = 4j(f,rc»7)» 

which has a free dependence on the variable 7, still possessing the initial condition 
70 = 0. In this case, the only invariance property is of the type (3.2) and there is no 
reduction other than that of second type. As above, one regards £ as a Lagrangian 
L{xiXi ^0)0) invariant under the isotropy group of 70 = 0. A simple concrete example 
of such a situation is when = (0,0,0) 6 M 3 and O = 50(3), which produces the 
framework for spin glass dynamics. Then, the momentum map associated to the residual 
50 (3)-symmetry (recall that 50(3) C ^(D, 50(3))) is J : T*F{V, 50(3)) so(3)*, 
J7"(a x ) = f v x -1 a! x At, and therefore yields (by Noether's theorem) the conservation law 

where the second equality follows from the general formula [Marsden & Ratiu(1994), 
formula (9.3.7)] 

|(Ad*a)=Ad^(d + ad* x _ 1 a) 

and from equations (3.8) with uq = (0, 0, 0) = n. As we have seen, in the case of liquid 
crystals, the two reductions are both possible (producing i\ and £2)- However, when i\ 
does not exist (e.g., for spin glasses), the only possible conservation law is 

V/fAd^uUo, (3.10) 



dt \J V V x Su , 

where i* is the dual of the Lie algebra inclusion i : o no °->- 0. In the case of spin glasses, 
n = (0,0,0) yields o no = so (3) = 0, so that i* reduces to the identity. On the other 
hand, for liquid crystals O no = 5 1 C 50(3) = O, so that the above conservation law is 
immediately implied by applying Neother's theorem to ii, as we did already in (3.4). 

The next section will apply this general setting to the case of microfluids. In this 
context the order parameter field is the molecule inertia tensor (microinertia) taking 
values in M = Sym(3), the space of 3 x 3 symmetric matrices. Both reductions above 
apply naturally in this context. 



12 



3.4 Reductions for micropolar media 



Micropolar media are continuum media in which the shape of each rigid particle may 
change in time, depending on the point in space. The molecule shape is given by an 
appropriate microinertia tensor, which also appears in the expression of the free energy, 
denoted by $. Then, the unreduced Lagrangian is given by 



p 



v 



Tr((z X 1 X) T X 1 x)^- / ®{X3oX \V(xjoX 



v 



where jo is the microinertia tensor and io := |Tr(j )^3 — jo or, equivalently, j = 
io — Tr(zo)/3. Upon repeating exactly the main steps as in the previous sections, one 
considers j G TiT), Sym(3)) as the order parameter field and, by defining j = x JoX~\ 
one obtains the first reduced Lagrangian 

■= I ^{xX~\xJoX~\^(xJoX~ 1 ))f^ = l I ®U,Vj)n, 

Jv 1 Jv Jv 

where 

Jo i-» XJoX^ 1 

defines the action of J r (V,SO(3)) on 3F(T>, Sym(3)). On the other hand, if j is constant 
in space, then we obtain the second reduced Lagrangian 



4(f,j,7) := / -&{XX \xjoX 1 ,[(Vx)x \xjoX 
Jv 

* Jv Jv 

where 7 := — (Vx)x _1 (with 70 = 0) and the Lie bracket [•, •] is the ordinary matrix com- 
mutator. Then, each of the above reduced Lagrangians produces the following equivalent 
sets of equations, respectively: 



d Si 2 

= 1/ x — 

at Su Su 



d Sit Sit 

— — = U X — 

at Su ou 
d t j + [j,u} = 0. 

Si 2 



S£i 



■ j 



(3.11) 



Sh 



. 1 . Sl.^f Si 2 

+ div h Tr 7 x — 

07 \ 07 



d tJ + [j,u] = 
{ d a + [7, u] + Vi> = , 7 = 0, 



(3.12) 



with the notation A» = e^A^. The explicit versions of the above equations are pre- 
sented in [Gay-Balmaz & Ratiu(2009)], where all the functional derivatives are also com- 
puted in detail. 
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While (3.11) are the well known equations for micropolar media [Eringen(1997)], the 
second set of equations (3.12) provide an alternative formulation of the same dynamics, 
as long as 7 = 0. When the latter initial condition is dropped, the two systems are not 
equivalent and (3.12) yield Eringen's formulation of micropolar liquid crystals, which 
accounts for disclination dynamics through the disclination density B in (1.1). Notice 
that the latter quantity possesses a coordinate-free definition in terms of the Yang- Mills 
curvature two-form 

B(u x , v x ) := d7(u x ., v x ) + 7(1^) x 7^) u x , v x E T X V , 

where d denotes the exterior differential. Then, if 70 7^ 0, the unreduced Lagrangian is 

£(X, X) = \ [ Tr {(toX^xfx-'x) [ ^(XJoX~\ X^X~ l + X7oX~>, (3.13) 
* Jv Jv 

with corresponding reduced expression 

e 2 (y,j,*y) = \( (3v)-v\i- [ *(j,7>- (3-14) 
* Jv Jv 

The discussion in the next section shows how this micropolar formulation recovers 
Ericksen-Leslie nematodynamics if the free energy \& equals the Frank energy (2.2). 

4 Comparing Eringen and Ericksen-Leslie theories 

In the previous section, Eringen's micropolar theory was shown to account for disclina- 
tion dynamics when 70 7^ 0. Now we shall show how Eringen's formulation of micropolar 
liquid crystals recovers EL nematodynamics, upon assuming that all molecules are uni- 
axial. In turn, this last constraint enforces the dynamics to neglect the presence of 
defects (i.e. 70 = 0), which otherwise would induce variations in the molecule shape. 

4.1 Micropolar theory for uniaxial nematics 

The assumption of uniaxial molecules enforces a microinertia tensor of the form 

j = J (I — n <g) n), 

which corresponds to i := |Tr(j)I — j = Jn<g>n. Here we assume that ||n|| 2 = 1. Then, 
this relation transforms Eringen's Lagrangian ^(f? j, 7), given in (3.14), to 

4(^n, 7 ) :=£ 2 (v, J(I-n®n),7) =- / |^xn|V- / ^( J(I - n ® n), 7 )/i, (4.1) 

^ Jv Jv 

thereby producing the equations (2.11), where the initial condition 70 = has not yet 
been imposed. Since we already know how equations (2.11), with Lagrangian (2.8), are 
related to EL nematodynamics, we need to choose a free energy \1/ which equals the 
Frank energy (2.2). Thus, we need to prove that there exists \l/ such that 

*(j, 7) = J{1 - n <g> n), 7) = F(n, n x 7) = F(n, Vn), (4.2) 
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where F is the Frank energy and the last equality follows from the relation Vn = n x 7, 
which is preserved by the dynamics (2.11). Notice that imposing Vn = n x 7 amounts 
to considering a subsystem of (2.11), unless 70 = 0. The next section shows how an 
appropriate free energy \1/ can be derived. 

4.2 The free energy 

This section shows how all terms in the Frank energy (2.2) can be rewritten in terms of 
the variables j = J (I — n <g> n) and 7, the latter being introduced through the invariant 
relation Vn = n x 7. Thus, the explicit expression for the micropolar free energy ^f(j, 7) 
of nematic media will be written after computing the micropolar expression for each term 
in the Frank energy (2.2). Some equalities are shown in detail in Appendix A. 

Twist. Using n ® n = I — j / J, we have 

n • V x n = -n • 7(11) + ||n|| 2 Tr(j) = — Tr (j'7) = — Tr (j'7 5 ) . 

J J 

where 7(11) = 7 ia n a , with a being the so (3) ~ M 3 -index, and *y s denotes the skew part 
of 7, i.e., 7 5 = (7 — 7 T ) /2, where we see 7 as a 3 x 3 matrix with components 7 io . 

Splay. We introduce the vector 7 = e abc 7 ac , defined by the condition 7-11 = Tr(ux7), 
for all u G IR 3 , where u x 7 is the matrix with components (u x 7)^ = (u x 7j) . We 
compute 

(div n) 2 = (7 • n) (7 • n) = 7 • (n ® 11)7 

= 7" Q--j/J)7 = Hf - j7" 31 

= 2 (Tr(j)/J - 1) Tr (( 7 A ) 2 ) - j Tr (j( 7 A ) 2 ) , 

where 7 A denotes the skew part of 7, i.e. "y A = (7 — 7 T ) /2 and where we used the 
equality 7 = — 27^. The latter can be shown by noting that we have the equalities 
Tr(7u) = -27 • u = -2 Tr(7u) for all u G M 3 . 

Bend. For all u G M 3 , we have 

(n x (V x n)) ■ u = — V n n ■ u = — (n x 7(11)) ■ u = — (u x n) ■ 7(11) 
= — un • 7(n) = — Tr((un) T 7n) 
= Tr ((n ® n)u 7 ) = Tr ((I - j/J)uy) 

= Tr(ux 7 )-iTr(ux ( 7J )) 
1 -> 

= 7 • u - -jJJ ■ u, 

so we get 

n x (V x n) = 7 — --yj 
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and therefore 



|n x (V x n) 



J 



13 ~1 



-2Tr 



{lj) A ~ 1 A 



Thus the micropolar form of the Frank free energy (2.2) is given by 



*0\7) = ^Tr( J7 ) + ^i( T r((7 A ) 2 ) (Tr(j) - J) - 2Tr(j(7 A ) 2 ) 



J 



J 



1 -^22 m_2 -^33 

2 J 2 



Tr 2 (j 7 )-^f Tr((( 7J ) A -J 7 A ) 2 ), (4.3) 



where we notice that the microinertia constant still appears explicitly. When the above 
free energy 7) is inserted in the micropolar Lagrangian £ 2 (u,j,^f), the equations 
(3.12) are completely equivalent to the Ericksen-Leslie equation (2.6). This is shown in 
the next subsection. 



4.3 Recovering Ericksen-Leslie nemato dynamics 

The relations j = J (I — n ® n) and Vn = n x 7 implied that the Frank energy has a 
micropolar formulation \I/(j, 7) = F(n, n x 7), to be used in the Lagrangian £ 2 (i/,j, 7). 
Here we restrict to the case 7 = 0, so that the equations (2.11) arising from £' 2 can now 
be transformed into a system of the type (2.3). In this way, the reduced Lagrangian 
iiiy, n, 7) transforms to a functional ii(u, n), by using F(n, n x 7) = F(n, Vn). Then, 
the resulting Lagrangian 

4(i/,n) = ^ f \uxnffi- f F(n,Vn)/x, (4.4) 



has to be inserted into the equations (2.3). Upon introducing the molecular field 



dn \dn ti J 
the variational derivatives are computed to be 



— — = — Jn x (n x u) = —Jv x [y x n) — h. 

ov on 



Equations (2.3) become 



J v x x (n x u)j — n x ^ J is x {y x n) + hj 



which reduce to 



J — (n x d t n) = J v x (n x 9 t n) — n x (Jvx d t n + h) . 



— J — ( n x (n x 1/ 
dtji + n x i/ = 
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The Jacobi identity produces thus the Ericksen-Leslie equations (2.6), thereby show- 
ing how Eringen's micropolar theory recovers Ericksen-Leslie nematodynamics. It is 
worth emphasizing that the solutions of the Ericksen-Leslie equations arising from the 
Lagrangian £[ lie on the following zero-level set 

5E 

-^•n = 0, 
ov 

thereby showing that there is no angular momentum Jnxh along the director field n. 
Again, this reflects the uniaxial nature of the nematic molecules described by Ericksen- 
Leslie theory. 



5 Two reductions for the Lhuillier-Rey theory 



The Lhuillier-Rey theory is an alternative description for liquid crystals of molecules with 
variable shape [Lhuiller & Rey(2004)]. In this description, the order parameter field is 
given by two components: the microinertia tensor j and the director field n. Because 
of the coexistence of these two order parameters, one can keep the expression of the 
Frank free energy (2.2) in the Lagrangian while allowing for a variable molecular shape, 
represented by the microinertia tensor appearing in the kinetic energy. The unreduced 
Lhuillier-Rey Lagrangian is 



£(XX)= / (x, X, X™o, XJoX \V(xno))/i 
Jv 

= o Tl ((ioX^x fx^X) V- / F(xn , V(xn ))// . 
^ Jv Jv 



(5.1) 



Again, this Lagrangian possesses two reductions, which are produced by the invariance 
properties 



£(XX)= / -&(x, X, X^o, XJoX 1 ,V(xn ))/i 
Jv 

= / ^(xX~\xno,XjoX~\ V(xn ))/i = £i(v, j, n) 
Jv 



(5.2) 



and 



£(x,x)= / ^(x,x,x n o,xjox l )v 

Jv 
Jv 



(5.3) 



where the second invariance arises naturally when no is constant in space, since one has 
Vn = (Vx)x _1 n. 

The first invariance yields the Euler-Poincare form of the Lhuiller-Rey equations 



v x 



d_5li 

dt Su 8u 
d t j + l], *>] = 
d t n + n x v = 0. 



-ft, 3 



8i x 
+ n x — 
on 



(5.4) 
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Reduction for the second group action yields the equivalent set of equations 




v 



(5.5) 



where we used the relation F(n, dn) = F(n, n x 7) (see appendix A). 

Notice that, the relations (Vj, Vn) = ([j, 7] , n x 7) determine an invariant sub- 
system, regardless of the initial conditions (j ,n ,7 ). Indeed, the possibility of an 
inhomogeneous initial condition 7 7^ extends the Lhuiller-Rey theory to account for 
disclination dynamics. Upon taking O = SO(3) and M = Sym(3) x S 2 , it is easy to see 
that the two Lhuiller-Rey formulations (5.4) and (5.5) follow directly from applying the 
general theory of Section 3 to the Lagrangian (5.1). 

It is interesting to notice that similar arguments to those in Section 4.3 show im- 
mediately how Lhuiller-Rey theory recovers Ericksen-Leslie nematodynamics. Indeed, 
while the relation Vn = n x 7 can be used to transform £ 2 into £1, the initial condition 



readily produces the Lagrangian (4.4), as was shown in Section 4.3, to recover the 
Ericksen-Leslie equations. 

6 Flowing liquid crystals 

This section extends the previous discussions to the case of flowing liquid crystals. In 
this case, upon denoting the fluid flow by the diffeomorphism rj £ Diff(I>), the gradient 
of the director field n = (x n o) ^ 1 — V*{x n o) is computed to be 

Vn = 77*(Vxn ) = ^((Vxx -1 ) X^o) = (v^^XX^ 1 )) (v*(x^o)) =: n x 7, 
where we have used standard properties of the push-forward 77* and we have defined 



A similar argument actually holds for any order parameter space M that is acted upon 
by a Lie group O. In this more general case, the previous relation reads as (3.5). This 
suggests that all the considerations in Section 3 have a natural correspondent in the case 
of liquid crystal flows. 

6.1 Euler-Poincare and affine Euler-Poincare reductions 

Upon restricting to incompressible fluid flows for convenience, the general form of the 
Lagrangian for flowing uniaxial liquid crystals is 



jo = J(I - n <g> n ) 



7 = -v*(^xx x ) • 
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so that £ is a functional of the type 

C : T(Diff vo i(£>) ® JT(V, O)) ->• K 

where the semidirect product structure reflects the fact that the microstructure vari- 
able x £ J-"(T>,0) is pulled around the domain V by the fluid flow r\ G Diff voi (X>) 
[Gay-Balmaz & Ratiu(2009)]. More precisely, the Lagrangian of a flowing uniaxial liq- 
uid crystal is 

£n (v,v,x,x) = o / \\v\\ 2 v + t; j [ IxlV- / ^((X^o) °^ _1 ,V((x^o) o^" 1 ))^, 

where F is the Frank free energy. Again, the fixed order parameter n Q appears in 
the Lagrangian density Jz? through both x n o an d its differential V(x^o)- Then, the 
Lagrangian L possesses the following invariance properties: 



v 



v 



V 



l 

Sf (?) o 7]' 1 , (xx' 1 ) ° J7 -1 , (xno) ° »7 -1 , V((xn ) o r/" 1 )) /i 
Jz? (u, n, dn) // = ^i(u, v, n) (6.1) 
and, if no is constant in space, 

£(x,x)=/ ■&( r liViX,X,XP>o, v(x%))^ 

-^('7 ° »7~\ (XX~ X ) ° »T\ (X^o) o r/" 1 , (?/•(( Vx) X _1 ))((x^o) ° J?" 1 )) /" 
J5f (u, i/, n, -7a/H) A* = ^( u > i/, ra, 7) , (6-2) 



v 



V 



V 

where one defines u := rjor/ -1 G X vo i(V), v := (XX _1 ) 07 7 _1 £ ^ r (^ > ) 0), n := (x ra o) or l 1 G 
F(P, M), and 7 := -r?*((V X ) x" 1 ) e fi 1 ^, 0). 

We notice that all the symmetry breaking arguments from Section 3 transfer to this 
case without essential changes. However, in the present case all variables are acted upon 
by the diffeomorphism rj G Diff vo i(X>), which affects the way variations are taken in the 
variational principles 

8 e 1 (u,u,n)dt = 0, 6 £ 2 (u, v, n, 7) dt = . 

J to J to 

Indeed, upon making use of the Lie derivative notation £, the variations 

5u = <9(\v + £ u w 

5u = d t oj — £^,v + £ u u + [lj, u] 

5n = —£ w n + u M (n) 

^7 = — -£w7 — V 7 a; 
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lead to the replacement all time derivatives in the equations (3.3) and (3.8) for motionless 
media by appropriate material derivatives, as is usual in the Euler-Poincare theory of 
fluid flows [Holm, Marsden & Ratiu(f998)]. Then, the resulting equations are 



D_6ii 
Dt~5^ 



— , Vv 
ov 



D5i x 8h 



5h 

— ■ dn 

on 



dp 



'SeA Dn 
' Sn r Dt 



(6.3) 



and 



p_5l2 
Dt 5u 



ov 



S£ 2 , S£ 2 
— ■ an + — 7 

on 07 



dp 



DS£ 2 S£ 2 J5£ 2 \^ A . J5£ 2 \ f Dn D>y\ 

Di^ + ad ^ = J [-^) +dw \Jt)> {-Dt'Dt) 



>mW,-VV) 



(6.4) 



where the material time derivative Dj Dt is expressed in terms of the Lie derivative as 

Dt dt + u ' 
while the diamond operator is defined by 



(/c 07,11) := - (k, £ v fl) . 
Also, Noether's theorems (3.4) and (3.10) become 



Dt V x Sv 







and — - % 
Dt 



v 



0. 



Remark 6.1 (Prom incompressible to compressible flows) Although in this sec- 
tion we focussed on incompressible fluid flows, extending to the compressible case re- 
quires no additional argument than those already present in [Gay-Balmaz & Ratiu(2009)]. 
Then, the mass density appears as an additional parameter in the Lagrangian C, whose 
symmetry properties (3.1) and (3.2) are accompanied by the relation p = V*P f° r the 
mass density p. This is also explained in [Holm, Marsden & Ratiu(1998)]. The dynam- 
ics of compressible liquid crystal dynamics can thus be obtained by simply combining 
the above results with standard principles in Euler-Poincare theory. 



6.2 Micropolar liquid crystals and Ericksen-Leslie theory 

This section extends the results of §4.3 to account for fluid motion. In the case of mi- 
cropolar liquid crystals, the order parameter field is j G Sym(3) and the Euler-Poincare 
Lagrangians arising from the two different reductions read 

£i(u,u,j) =i / |u|V+o / (» / $(j,dj>, 

^ Jv £ Jv Jv 
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and 

where compatibility of the two reductions enforces the energy relations 

$(j,dj) = $(j, [7,7]) = *(j,7)- 

Then, upon specializing to uniaxial nematics, the microinertia tensor is given by 

j = J(n (g) n — I) , 

and the corresponding free energy must transform as 

7) = \P(J(n <g> n — I), 7) = F(n, n x 7) = F(n, Vn) . 

The expression for the free energy ty(J, 7) is given in (4.3). In conclusion, the Lagrangian 
£2 is transformed to 

4(u,i/,n) = - / \u\ 2 /i+-J \vxn\ 2 fi- / F(n,Vn)/i, 

which in turn produces the equations (6.3) in the form 

— -r 1 = -Vi/ • -r 1 - Vn ■ dp 

Dt ou ov on 

-D < 54 < Dn 

— v x — = n x — , — — = v x n. 

oV 8u 5n Dt 

Upon inserting the variational derivatives 

— i = u, — - = — J n x (n x v) , — - = — h — Jv x [y x n) 
ou ov on 

and by repeating analogous steps to those in §4.3, one obtains the hydrodynamic Ericksen- 
Leslie equations 

f / 

— + (u • V)u = -Vp - d l ( Vn 



<9t \ cfrij 

D 2 n / , £> 2 n\ 
J— -2 n.h + in.- )n + h = 0. 



Again, this proves that Eringen's micropolar theory of liquid crystal flows recovers 
Ericksen-Leslie theory of uniaxial liquid crystal flows. 

A Expression of the Frank energy 

We show that the relation Vn = —7 x n allows the Frank free energy F(n, Vn) to be 
expressed only in terms of 7(x) and n(x), so that F(n(x), Vn(i)) = \l/(n(x), 7(2?)). 
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Twist. First, we compute in coordinates: 

n • V x n = n l e ijk djn k = -n l e ijk e kab j a n b = -n l n b (e kij e kab ) 7^ 
= - n l n b (5 ia S jb - S ib S ja ) 7^ = -n a nrf a + n b n b ^ a a . 

Coordinate free computation: 

n ■ (V x n) = — Tr (fi2(Vn) A ) = - Tr (nVn) 

= — Tr (n(n x 7)) = — Tr (nn7) 
= — n ■ 7(n) + ||n|| 2 Tr(7). 

Splay. Coordinate computation: 

(divn) 2 = {diUi) 2 = (e iab Y a n b ) 2 = e iah f a n b e jcd j J c n d 
= (eiafe7a e jcd7c) n dn b 

Coordinate free computation: We define the vector 7 by the equality 
7-u = Tr(ux 7 ), i.e., 7 = e io6 7^ or 7 = -27^ 

and we compute 

(divn) 2 = Tr(Vn) 2 = Tr(n x 7) 2 = (7 • n)(7 • n). 

Bend. Coordinate computation: 

n x (V x n) = e ijk n J e kab d a n b = - (e kij e kab ) n J (e bcd ^n d ) 
= - (SiJjb - 5 ib 5 ja ) n ] (e bcd ^n d ) 
= ~ Siadjhn 1 (e bcd j^n d ) - 5 ib S ja n 3 (e bcd ^n d ) 
= - (tjafll + ticdll) n d n 3 . 

Coordinate free computation: For all u 6 M 3 , we have 

(n x (V x n)) • u = — V n n ■ u = — (n x 7(n)) • u = — (u x n) • 7(n) 
= — un • 7(n) = — Tr((un) T 7(n)) 
= Tr ((n (8) n)u7) , 

so we get 



n x (V x n) = 7(n ® n) = — n x 7(n) 

and therefore, 

||n x (V x n)|| 2 = ||n x 7 (n)|| 2 = ||n|| 2 ||7(n)|| 2 - (n ■ 7 (n)) 2 
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Putting all these results together, we conclude that there exists a function \I/(n(x), 7(x)), 
such that \l/(n(x), 7(x)) = F(n(x),Vn(i)), that is, 

(n(x), 7 (x)) = -K 2 (||n|| 2 Tr( 7 ) - n • 7 (n)) + ^11(7 • n) 2 

+ i/i 22 (||n|| 2 Tr(7)-n-7(n)) 2 
+ ^3 3 (||n|| 2 ||7(n)|| 2 -(n-7(n)) 2 ). 

References 

[Chandrasckar(1992)] Chandrasckhar, S. [1992], Liquid Crystals, Cambridge University Press, Cam- 
bridge. 

[dcGcnnes(1971)] de Gennes, P. G. [1971], Short range order effects in the isotropic phase of nematics 
and cholesterics. Mol. Cryst. Liq. Cryst., 12, 193-214. 

[Dzyaloshinskii & Volovick(1980)] Dzyaloshinskii, I. E. and Volovick, G.E. [1980], Poisson brackets in 
condensed matter systems, Ann. Phys. 125, 67-97. 

[Ellis et al.(2010)] Ellis, D. C. P., Gay-Balmaz, F., Holm, D. D., Putkaradze, V., and Ratiu, T. S. 
[2010], Symmetry reduced dynamics of charged molecular strands, Arch. Rational Mech. Anal. 
197(3), 811-902. 

[Eringcn(1997)] Eringen, A. C. [1997], A unified continuum theory of liquid crystals, ARI - An Inter- 
national Journal for Physical and Engineering Sciences 50, 73-84. 

[Gay-Balmaz, Holm & Ratiu(2011)] Gay-Balmaz, F., Holm, D. D., and Ratiu, T. S. [2009], Variational 
principles for spin systems and the Kirchhoff rod, J. Geom. Mech., 1(4), 417-444. 

[Gay-Balmaz & Ratiu(2009)] Gay-Balmaz, F. and Ratiu, T. S. [2009], The geometric structure of com- 
plex fluids, Adv. Appl. Math., 42(2), 176-275. 

[Gay-Balmaz & Tronci(2010)] Gay-Balmaz, F. and Tronci, C. [2010], Reduction theory for symmetry 
breaking with applications to nematic systems, Physica D, 239(20-22), 1929-1947. 

[Holm(2002)] Holm, D. D. [2002], Eulcr-Poincare dynamics of perfect complex fluids, in Geometry, 
Dynamics and Mechanics: 60th Birthday Volume for J.E. Marsden. P. Holmes, P. Newton, and 
A. Wcinstcin, eds., Springer- Verlag. 

[Holm & Kupcrshmidt(1988)] Holm, D. D. and Kupcrshmidt, B. A. [1988], The analogy between spin 
glasses and Yang-Mills fluids, J. Math. Phys. 29, 21-30. 

[Holm, Marsden & Ratiu(1998)] Holm D. D., Marsden, J. E., and Ratiu, T. S. [1998], The Eulcr- 
Poincare equations and semidirect products with applications to continuum theories, Adv. in 
Math., 137, 1-81. 

[Lcslic(1979)] Leslie, F. M. [1979], Theory of flow phenomena in liquid crystals, in Advances in Liquid 
Crystals, vol. 4, (G. H. Brown, cd.) Academic Press, New York pp. 181. 

[Lhuiller & Rey(2004)] Lhuillicr, D. Rey, A. D. [2004], Nematic liquid crystals and ordered micropolar 
fluids, J. Non-Newtonian Fluid Mech. 120, 169-174. 

[Marsden & Ratiu(1994)] Marsden, J. E. and Ratiu, T. S. [1994] Introduction to Mechanics and Sym- 
metry. Springer- Verlag. 

[Rey & Denn(2002)] Rey, A. D. and Demi, M. M. [2002] Dynamical phenomena in liquid-crystalline 
materials. Ann. Rev. Fluid Mech. 34, 233-266. 



23 



